Calculation of Free-Atom Fractions in Hydrocarbon-Fueled Rocket 

Engine Plume 


Abstract 

Free atom fractions ((3) of nine elements are calculated in the exhaust plume of CH 4 - 
oxygen and RP-1 -oxygen fueled rocket engines using free energy minimization method. 
The Chemical Equilibrium and Applications (CEA) computer program developed by the 
Glenn Research Center, NASA is used for this purpose. Data on variation of [3 in both 
fuels as a function of temperature (1600 K - 3 100 K) and oxygen to fuel ratios ( 1 .75 to 
2.25 by weight) is presented in both tabular and graphical fonns. Recommendation is 
made for the P value for a tenth element, Palladium. The CEA computer code was also 
run to compare with experimentally determined P values reported in literature for some of 
these elements. A reasonable agreement - within a factor of three - between the 
calculated and reported values is observed. Values reported in this work will be used as a 
first approximation for pilot rocket engine testing studies at the Stennis Space Center for 
at least six elements Al, Ca, Cr, Cu, Fe and Ni - until experimental values are generated. 
The current estimates will be improved when more complete thermodynamic data on the 
remaining four elements Ag, Co, Mn and Pd are added to the database. A critique of the 
CEA code is also included. 

Introduction 

Development and testing for flightworthiness of a rocket engine is critically dependent 
upon monitoring the exhaust plume. Ten elements (Ag, Al, Ca, Co, Cr, Cu, Fe, Mn, Ni 
and Pd) are spectroscopically peculiar to the materials used in rocket engines 1 and are of 
special interest to the plume diagnostic team at the John C. Stennis Space Center, 
Mississippi. Detection, identification and quantification of these elements is a part of 
rocket engine health monitoring program which can indicate the rate of progressive 
malfunction of specific components during the engine testing. Proper actions may then 
be initiated to avoid a disaster. 

Real time detection of these elements is accomplished by monitoring and analyzing the 
spectral signatures from the atoms of these elements in the exhaust plume of a rocket 
engine. The total rate of discharge of each element in the plume is calculated by 
appropriately applying a quantity called P for each element. P is defined as the mole 
fraction of an element that is in the atomic form versus all its forms in the flame e.g. 
molecular, ionized, oxide, hydride, hydroxide etc. The spectroscopic data is regressed 
through several layers of computer codes to quantify the rate of wear and identify the 
particular engine components as the source of the discharge of a particular element. 

Details of the methodology of rocket engine health monitoring using spectral analysis, 
computer codes, scale model tests and ground tests for flightworthiness of the space 
shuttle main engines have been extensively published by Tejwani et.al . 2, 3 ' 4 ’ 5 


RELEASED - Printed documents may be obsolete; validate prior to use. 



2 


NASA has been investigating the suitability of high energy density hydrocarbon fuels for 
the new generation reusable space shuttle engines. Two main hydrocarbons fuel 
candidates are CH4 and RP-1. A comprehensive literature 6 search was conducted in 
summer 2005 and a compilation was made of the P values for the above mentioned ten 
elements in flame environments similar to the hydrocarbon-oxygen propelled rocket- 
engines. The environments reported in the literature were acetylene-oxygen, acetylene- 
nitrous oxide and hydrogen-nitrous oxide flames. P values generated from the current 
space shuttle engine tests using liquid hydrogen-liquid oxygen fuel system were also 
included in the above study for comparison purposes. 

We concluded in that study that the reported P values could not be used for CH4 or RP-1 
flame conditions, because of the sensitivity of P values to flame composition, uncertainty 
in the details of experimental conditions and the large variations in the P values reported 
by different investigators. Alternative experimental and theoretical methods were 
recommended to generate P values for the CH4 and RP-1 flames at oxygen / hydrocarbon 
ratios of interest to us. It was also recommended that because our primary interest is in 
the equilibrium composition, a thermodynamic approach rather than a kinetic approach 
should be used. 

Kinetic Approach vs. Free Energy Minimization Approach 

Equilibrium composition of a reaction mixture may be calculated by both kinetic and 
thennodynamic methods. The kinetic approach utilizes the equilibrium constants of the 
various elementary reactions occurring in the mixture at equilibrium: 


aiA + biB <-> CiC + diD 

Ki=([C ] c1 [D] dl ) / ([A] al [B] bl ) 

( 1 ) 

a 2 A + eiE <-> fF + giG 

K 2 = ([F] fl [G] gl )/([A] a 2 [E] el ) 

( 2 ) 

a 3 A + f 2 F hiH + liL 

K 3 = ([H] hl [L] n )/([A ] a3 [F f) 

( 3 ) 

b 2 B + e 2 E <— > j 1 J + 1 2 L 

K 4 = ([J] jl [L] 12 )/([B ] b2 [E] e2 ) 

( 4 ) 


K n = (n) 

Where A, B, C etc. are the species generated from the starting compounds in the reaction 
mixture, al, a2, bl, b2, etc. are their coefficients in the reactions, [A], [B], [C] etc. are 
their molar concentrations in the equilibrium mixture and Ki, K 2 , K 3 etc. are the 
equilibrium constants for the reactions. Thus, a set of equations accounting for all 
possible interactions of all possible species may be considered, each with its own 
equilibrium constant at the appropriate reaction temperatures. The equilibrium 
concentrations of m species can be obtained by simultaneous solution of n equations in m 
species (where m < n). 
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By comparison, free energy minimization (FEM) method utilizes the fact that the Gibbs 
free energy of a system is at a minimum when it is at equilibrium. The change in Gibbs 
free energy of a component is defined by: 

Gi = H; - T * S; 

Where, 

Gi = Gibbs free energy of component i. 

H; = Enthalpy of component i. 

S; = Entropy of component i. 

The enthalpy and entropy of a component are functions of temperature and are defined in 
tenns of specific heat at constant pressure C p ° by (component subscript i is not used to 
reduce clutter): 

k 

C p °/R = I a k *T' k 
H°/RT = | C p ° dT/ RT 
S°/R = J (C p ° / RT) dT 

R is the Universal Gas Constant. Thus, if the temperature dependence of C p ° for a 
component can be described as a polynomial with seven tenns (k = 7), the enthalpy and 
entropy can then be defined as follows: 

C p °/R = ai*T " + a 2 *T 1 +a 3 + a 4 *T + a 5 *T" + a6*T 3 + a 7 *T 4 

H°/RT = -ai*T' 2 + a 2 *T' 1 *lnT + a 3 + a 4 *T/2 + a 5 *T 2 /3 + a 6 *T 3 /4 + a 7 *T 4 /5 + bfiT 

S°/R = - ai*T 2 /2 - a 2 *T 1 + a 3 *lnT + a 4 *T+ a5*T”/2 + a6*T 3 /3 + a 7 *T 4 /4 + b 2 

bi and b 2 are integration constants. Thennodynamic data on only those chemical species 
existing at equilibrium is needed to calculate the equilibrium composition. As a result, the 
need to know the exact composition of the initial fuel and the exact reactions between the 
numerous species is circumvented. 

The free energy of a system or the change in free energy of a system starting from a 
standard state, which contains n species at molar concentration c, at equilibrium, can be 
calculated from the individual AG; thus: 

n 

AG = I Ci AG, 


It is Ci, the molar concentration of species i at equilibrium, is the quantity we strive to 
estimate. 
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It should be pointed out here, that the kinetic and thermodynamic approaches should and 
must yield the same equilibrium composition provided the data base is complete. The 
equilibrium constant K and the change in the free energy AG of a system are related by 
the well known equation: 

AG = -RT In K 

Thus if one is obtained from experimental or theoretical means, the other may be 
calculated. 

To recapitulate, if all reactions occurring in the flame and their equilibrium constants are 
known, one can calculate the equilibrium concentration of species in the flame. For 
simple flames this approach is suitable, however, for complicated systems, this approach 
is not efficient because all flame reactions may not be sufficiently characterized in terms 
of the species involved, and activation energies and temperature dependent reaction rate 
coefficients for all these reactions may not be available. Consequently, most investigators 
limit the consideration to the interaction of a handful of species of highest concentrations. 
For example, a kinetic model for combustion of acetylene considers thirty two species 
and one hundred and six reactions ! Another review article on combustion of various 
hydrocarbons reports fifty three species, three hundred twenty five elementary reactions 
and associated rate coefficient expressions and thermochemical parameters for 
combustion of CFG and twelve hundred reactions and one hundred twenty seven species 
for ethane and propane combustions. Given the computing power available now-a-days, 
numbers of species and reactions may no longer be an issue, however, the problem 
becomes significantly more complex in the case of mixed fuels such as kerosene, jet fuels 
such as RP-1 and rocket propellants where the fuel compositions is variable. The 
presence of non-ideal species and condensed phase introduces further difficulty in the 
kinetic approach. By contrast, in the FEM method, thennodynamic data on only those 
species existing at equilibrium is needed to calculate the equilibrium composition. As a 
result the need to know the exact composition of the initial fuel and the exact reactions 
between the species is circumvented. A surrogate formula may be used for mixed fuels in 
the FEM approach and the presence of ionized and condensed species at equilibrium 
conditions does not pose any problem. 

The free energy minimization method has been used in the past for combustion 
reactions 9 ' 13 . Some authors have used the classical thermodynamic approach 9 ' l0 ' 11 and 
some others have suggested a statistical mechanical approach ’ . In this report we have 
used a computerized method known as the Chemical Equilibrium and Application (CEA). 
This FORTRAN based computer code was developed and refined by Gordon and 
McBride at the Glenn Research Center NASA 14 ’ lx 16 over a number of years and has 
been used extensively since then. The CEA program requires the following information 
on the species: heat of formation at 298. 15 K, phase (gas / condensed) of the specie, its 
molecular weight, temperature coefficients (ai through a 7 ) in the equation for heat 
capacity at constant pressure C p ° and the temperature intervals over which they are valid, 
assigned enthalpy H°(T) at a reference temperature T to correct for non-ideality, and 
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integration constants (bi, b 2 ) in the expressions of enthalpy and entropy (expressed as 
dimensionless quantities H°/RT and S°/R respectively). 

A brief comment on whether equilibrium is actually achieved in the flame reactions is as 
follows. The presumption of thennodynamic equilibrium implies that the rate of reaction 
between various species is fast compared to the residence time of these species in the 
flame. Attainment of equilibrium requires that the particles in the system have undergone 
a sufficient number of collisions with each other. According to Alkemade 17 , calculations 
show that the translational and rotational energies of particles equilibrate very quickly - 
within 10 ‘ 8 sec at atmospheric pressure. Although, vibrational energy equilibration is 
slower in comparison with translational and rotational energies, any serious lag in 
equipartition of energy is not expected in this time frame. It is, therefore, safe to assume 
that local thermodynamic equilibrium is achieved in the flame reactions. 

An interesting feature of the CEA code, and indeed equilibrium compositions at specified 
temperature and pressure should be pointed out here. If the oxygen/carbon and 
oxygen/hydrogen molar ratio in the feed are kept constant, the equilibrium distribution of 
species does not change whether the starting compound is C 2 H 4 or C 3 H 6 or C 4 H 8 and so 
on - because these compounds also have a constant hydrogen/carbon molar ratio. Stated 
in general terms, all reactions characterized as 

C n H m + p O 2 = various C, H, O species 

will yield the same equilibrium composition for all feed compositions that have the same 
n/p and m/p ratios. The coefficients n, m and p may be fractions. 

The two determining factors for making this statement are: the relative thennal stability 
of species at the specified temperature and pressure and the mass balance on elements. 
Different heats of combustion for different reactants does not affect the equilibrium 
composition because temperature and pressure are specified. This is also the basis for 
using a surrogate formula for RP-1 in the absence of a “true” formula. 

RP-1 surrogate formula 

The hydrogen to carbon (H/C) ratio for RP-1 reported in literature has varied over a wide 
range ’ ’ . Latest RP-1 specification indicates a minimum hydrogen content of 13.8 
% by mass. This would correspond to an H/C ratio of 1.91. RP-1 in the CEA data base 
has the elemental formula CH1.95. Wang used a surrogate formula of C12H24 for RP-1 
(C/H ratio of 2.0) for fluid dynamic modeling of the combustion process. We expect that 
efforts to decrease the aromatic fraction, (H/C ratio =1) in RP-1 to avoid soot formation 
tendency 19 , and increase naphthalene (H/C ratio =2) and paraffin fractions (H/C ratio >2) 
will move the desired RP-1 formulation toward an H/C ratio of two. Using an H/C ratio 
of two or a surrogate formula of C 12 H 24 keeps our results in confonnity with the fluid 
modeling efforts. 
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A further point is noted here regarding the formulation of RP-1. Refinery processes such 
as catalytic cracking, alkylation etc, are very corrosive and it is likely that the RP-1 fuel 
may contain Fe, Cr in trace quantities. Also, Ni, Co, Ag are commonly used as catalysts 
in the oil refining process, therefore, these elements may also be present at parts per 
million level concentrations. Therefore, it is necessary to establish a base line 
concentration of these elements in RP-1. The CR-101 thruster tests 5 with RP-1 will 
provide an opportunity to make such detennination. The thruster will operate on pressure 
feed without the aid of any rocket engine turbomachinery. 

Case Studies Run 

Two cases reported in literature were chosen to run on CEA so that a comparison with the 
reported results can be made. Then CH 4 and RP-1 cases of interest to us were run using 
the CEA code. 

The combustion of CH4 and RP-1 is represented thus: 

CH 4 + 2 O 2 = CO? + 2 H 2 O O/F ratio = 4.0 for complete combustion (cp = 1) 

C 12 H 24 + 18 O 2 = 12 CO 2 +12 FEO O/F ratio = 3.43 for complete combustion (cp = 1) 

The choice of oxygen / fuel ratios (O/F ratios) for CH 4 , O 2 system was made to include 
the ratios of 2.0 and 2.5 used in initial studies on CH 4 thruster currently being conducted 

Table -1 Case Studies Run on the CEA Code 


Comparison with cases reported in literature 


O/F ratios by weight 
(Equivalence ratios) 

Temperatures K 

Pressure 

atm 

Elements # 

C 2 H 2 , 0 2 

1.368, 1.758 
(0.445,0.571) 

2750 

1.0 

Ca, Cr, Cu, Fe 

h 2 , o 2 

6.0** 

3097 

2.455 

Cr, Cu, Fe, Mn, Ni 

Cases in the current study 

CH 4 , 0 2 

2.0,2.25,2.5 
(0.5, 0.563, 0.625) 

1600, 2000, 2400, 
2600, 2800, 3000, 
31000 

2.0 

Ag, Al, Ca, Co, Cr, 
Cu, Fe, Mn, Ni 

RP-1,0 2 

1.75,2.0,2.25 

(0.51,0.583,0.656) 

1600, 2000, 2400, 
2600, 2800, 3000, 
31000 

2.0 

Ag, Al, Ca, Co, Cr, 
Cu, Fe, Mn, Ni 

RP-1, S*, 
0 2 

1.75,2.0,2.25 

(0.51,0.583,0.656) 

1600, 2000, 2400, 
2600, 2800, 3000, 
31000 

2.0 

Ag, Al, Ca, Co, Cr, 
Cu, Fe, Mn, Ni 


# All elements were present at 1 ppm by mole concentration 

* 30 ppm Sulfur by weight 

**Space shuttle main engine test conditions 
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at the Stennis Space Center. These correspond to equivalence ratios cp = 0.5 and 0.625 
respectively. The O/F ratios for RP -1 were chosen arbitrarily in the vicinity. 

The CEA program was run as a constant temperature, constant pressure (t, p) problem. Pd 
could not be included in these calculations because it is not in the CEA data base. In the 
case of CH4 and RP -1 feeds without sulfur, the mole fraction of each of nine elements in 
the hydrocarbon feed was 1 . 0 x 1 O' 6 , and that of the hydrocarbon fuel itself was 99 . 9991 . 

In the case of RP -1 fuel with sulfur, the mole fraction of sulfur was 0.0157 corresponding 
to 30 ppm by weight in the fuel of fonnula C12H24. Consequently, the mole fraction of the 
RP -1 fuel was 99 . 9834 . Oxygen was input as 100 % pure oxygen. The input compositions 
for various cases are presented in Table- 2 . 

Ions, electron gas and all possible species combinations of C, H and the elements in the 


Table - 2 Input File for CEA Program for Various Case Studies 


Case 

C 2 H 2 , 0 2 

h 2 , 0 2 

CH 4 , 0 2 

RP- 1 , 0 2 

RP- 1 , S, 0 2 

O/F ratio 

Variable, see r 

fable 1 

0 2 

100.0 

100.0 

100.0 

100.0 

100.0 

ch 4 

99.9996 

99.9995 

99.9991 



RP -1 




99.9991 

99.9834 

S 





.0157 

Ag 



.0001 

.0001 

.0001 

A 1 



.0001 

.0001 

.0001 

Ca 

.0001 


.0001 

.0001 

.0001 

Co 



.0001 

.0001 

.0001 

Cr 

.0001 

.0001 

.0001 

.0001 

.0001 

Cu 

.0001 

.0001 

.0001 

.0001 

.0001 

Fe 

.0001 

.0001 

.0001 

.0001 

.0001 

Mn 


.0001 

.0001 

.0001 

.0001 

Ni 


.0001 

.0001 

.0001 

.0001 


data base were included in the equilibrium mixture. The program output was in tenns of 
mole fraction of various species in the equilibrium mixture. The lower limit on the 
concentration of species reported in the output was reduced to lxlO' 15 mole fraction from 
the default value of lxl 0' 5 . [3 for each element was calculated using the following 
equation: 

[3 = (mole fraction in atomic form)/ E (mole fractions of the element in all forms) 

Air entrainment in the plume was not taken into consideration. Therefore, no nitrogen 
species were assumed to be present in the equilibrium mixture. The calculated (1 values 
for all cases are presented in Tables 4 through 7 and Figures 1 and 2 . 
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Discussion 

a. Critique of the CEA Program 

The CEA program, especially with its current graphical user interface, is an excellent tool 
for calculations of equilibrium composition. The program is essentially a mass balance on 
atoms of all elements in the feed distributed among the various species containing these 
elements in the data base, while seeking the minimum of the free energy for the whole 
system at the same time. Any computational program, including CEA, is completely 
dependent upon the quality and comprehensiveness of its data base. If all relevant forms 
of an elements e.g. ions, radicals and compounds with other elements, are not included in 
the data base, an erroneous mass fraction for free atoms of that element will result. 

Further error may be introduced if the temperature dependence of enthalpy and entropy 
for each species is not available for rigorous energy minimization calculations. 

Table 3 lists the species that are in the data base for each element of interest to us in the 
hydrocarbon combustion environment. Most species in the data base have temperature 
dependent data up to 6000 K, and some up to 20000 K. Data on remaining species is 
valid for temperatures less than 2500 K, they are crystals forms and are not expected to 
exist above their melting temperatures. 

The data in CEA is taken from the JANAF tables and other reputed sources, therefore, 
it is considered to be of good quality. Pd is absent from the database. Also as noted in a 
later section, some species of interest to us, such as Pd species (e.g. Pd+, Pd 2 H, PdFE, 
PdO), Ag species (AgOH, Ag 2 0, Ag 2 0 2 , Ag 2 0 3 , Ag 3 0 4 , Ag 3 0, Ag 2 S), Co species (CoO, 
Co 3 0 4 , and CoS0 4 ) and Ni species (NiH and Ni 2 0 3 ) are not in the data base. Similarly a 
limited number of Mn species are present in the data base. The extent, to which inclusion 
or absence of these species would affect the equilibrium distribution of an element in 
various forms at temperatures of interest to us, is unknown. There is a provision to add 
data on new compounds to the database. However, data on all possible species of interest 
to us is not yet available in the literature. The data base on Al, Ca, Cr, Cu, Fe and Ni 
species appear to be adequate; therefore, the results for these elements should be more 
readily acceptable. 

b. Comparison of CEA results with experimental values reported in literature 

In order to assess the validity of the results predicted by CEA program, two cases were 
chosen from the published literature. One of the reports is an excellent review and 
comparison of data in open literature on flame spectroscopy and of relevance to us 
because it is a hydrocarbon in oxygen flame (C 2 H 2 -0 2 flame diluted in He). It should be 
noted here that the objective of a majority of work reported in literature has been to 
artificially increase the free atom fractions in the flame by addition of K, Cs etc. as 
ionization suppressants. This practice reduces the number of cases available for valid 
comparison by a significant margin, because, only those cases where no foreign agent 
was added to the flame may be used. It is expected that in the example used, Helium does 
not participate in any of the flame reactions; therefore, P values for the elements remain 
unaffected and a fair comparison may be made between the calculated and 
experimentally determined values. 
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Table-3 List of Species Included in the Free Energy Minimization Calculations 



e- 


CH20H 

EE 

C4H8,1 -butene 

■Mil 


eeei 


ETiH 



Ag 

58 

CH20H+ 

mi 

C4H8,cis2-buten 

■Ml 


Kfli. 

CuS04(cr) 

KM 

Mn(d) 


Ag- 

59 

(CH3COOH)2 

HI? 

C4H8,cyclo- 

■Ml 

CaH 

t 

Fe 

KiM 

Mn(L) 

4 

Ag+ 

ee 

(HCOOH)2 


C4H8, isobutene 


CaH2(a) 


Fe- 

mi 

Ni 

m 


61 

CH3CHO,ethanal 


C4H8,tr2-butene 


CaH2(b) 


Fe+ 


Ni- 

6 

Ag(L) 


Mi— 




CaH2(L) 


Fe(a) 


Ni+ 

m 

AL 


CH3COOH 


C4H9,n-butyl 


CaO(cr) 


Fe(c) 

fWi 

Ni(cr) 

8 

AL2 

mi 

CH30 




CaO(L) 

mxj 

Fe(CO)5 


My 


AL- 

m 

CH302CH3 


C4H9,t-butyl 


CaOH 


1 1 llilil HIM 




AL+ 

66 

CH30CH3 


C5H10,1-pentene 


CaOH+ 


Fe(d) 


Ni3S2(b) 


ALO 

67 

CH30H 




CaS 

vm 

Fe(L) 


I3B3- 


AL(cr) 

68 

CH300H 




CaS(cr) 

ee 

Fe(OH)2 


liESSiHiiiil 


AL(L) 


CH4 




CaS(L) 




NiO 


AL(OH)2 

70 

COOH 




CaS04(l) 


iMiS! 


NiS 

HE 


El 

COS 

Eh 

C5H12,n-pentane 

183 

CaS04(ll) 

239 

Fe.9470(cr) 


NiS(a) 

m 



CS 

■Mil 

C5H6,1,3cyclo- 

■EEI 


EE 

Fe.9470(L) 

Kilil 

NiS(b) 

m 

AL2C2 

m 

CS2 

Bfefel 

C5H8,cyclo- 

■l;M 

CH 


Fe2(S04)3(cr) 

Kiri 

NiS(L) 

he 

AL20 

74 

C2- 

be 

C6H10,cyclo- 

■ Mil 

CH+ 


Fe203(cr) 

EEI 

NiS2(cr) 

HE 

AL20+ 

E 

C2H 

■Ml 

C6H12, 1-hexene 

■I.Vj 

CO 

ee 

Fe304(cr) 

KM 

NiS2(L) 

Til 

AL202 

mi 

C2H2, acetylene 

■IcVJ 

C6H12,cyclo- 

■Mil 

Co 

eq 

Fe304(L) 

ee 

0 

El 

AL202+ 

■MU 


■Mil 

C6H13,n-hexyl 


Co- 


FeO 


0- 


AL203 

78 

C2H3, vinyl 

134 

C6H14,n-hexane 


Co+ 


FeS(a) 


0+ 


AL203(a) 

79 

C2H4 


C6H2 


Co(a) 


FeS(b) 


02 

m 

AL203(L) 

eh 

C2H40,ethylen-o 


C6H5, phenyl 


Co(b) 


FeS(c) 


02+ 

EE 

AL2S 

81 

C2H5 


C6H50,phenoxy 


C°(L) 


FeS(L) 


OH 


AL2S2 

m 

C2H50H 


C6H50H, phenol 


CO+ 

ee 

FeS2(cr) 


OH+ 



m 

C2H6 


C6H6 


C02 


FeS04(cr) 




, Ej.i i'SF 

MM 

C20 




C02+ 

vm 

H2- 


02- 


• i-w ■ 

HE 

C2S2 




Cr 


H20 


03 


l | M i i| 

EE 
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Cr- 


H20(cr) 


OH- 


ALC 

m 

wikmaH 
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Cr+ 


H20(L) 




ALC2 

88 

C3H4,allene 
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C7H7, benzyl 

fcWil 

CrO 

ee 

H202 

ee 

S 


ALH 

89 

C3H4,cyclo- 

EE 

C7H8 
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Cr(cr) 

E h 

H2S 

ee 

S- 

HQ 

ALH2 

Mil 

C3H4,propyne 

■EM 

C7H80,cresol-mx 

WiW 

Cr(L) 

ee 

H 

eei 

S+ 

US 

ALH3 

91 

C3H5,allyl 

%Vrl 

C8H10,ethylbenz 

BE 

Cr203(l) 

ee 

H- 

ee 

S(a) 


ALH3(a) 

92 

C3H6,cyclo- 

■EM 

C8H16,1-octene 

Effil 

Cr203(l') 

ee 

H+ 

ee 

S(b) 

m 

ALO- 

KB 

C3H6, propylene 

■E!;' 

C8H17,n-octyl 

WiM 

Cr203(L) 

mi 

H2 

Eh 

S(L) 


ALO+ 

94 

C3H60, acetone 

150 

C8H 18, isooctane 

EE 

Cr02 


H2+ 


S2 


AL02 

B33 

C3H60,propanal 


C8H18,n-octane 


Cr03 


H20+ 


S2- 

EE 

AL02- 

96 

C3H60,propylox 




Cr03- 


H30+ 


S20 

ms 

ALOH 

EE 



C9H19,n-nonyl 


Cu 


HCO+ 


S3 

EE 

ALS 

98 

C3H7,n-propyl 


C10H21,n-decyl 


Cu- 


H2S04 


S4 

EE 

ALS2 


C3H8 


C10H8,naphthale 


Cu+ 

EE 

1 1 iHlIM 


S5 


C 


C3H80,1 propanol 




l. '.'il*. . ‘V- 


HALO 


S6 


C(gr) 


C3H80,2propanol 


MhBMHWWHlM 




HAL02 


S7 


c- 


C302 
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Ca 


mmmamm 


HCCO 


S8 

KQ 

c+ 


C30S 


Ca2 


Cu2 



fm 

SH 

EEI 

C2 

UiQ 

C3S2 

■rail 

Ca(a) 

EE 



HCO 

ee 

SH- 

EH 

C2+ 

■MS 

C4H 10, isobutane 

■Ml 

Ca(b) 

Eb 

Cu20(L) 

ee 

HCOOH 

ee 

SO 

wail 

C3 

■Mil 

C4H10,n-butane 

■r*i 

Ca(L) 

be 

Cu2S(a) 

Vrl 

H0(C0)20H 

ee 

SO- 

mi 

C4 

■ItW 


■r*i 

Ca+ 

ee 

Cu2S(b) 


H02 

Eil 

S02 

H 

C5 

■Mil 

C4H4,1,3-cyclo- 

■Ml 

CaO 

fcMil 

Cu2S(c) 

fiffl 

H02- 

ee 

S02- 

wan 

CH2 

IfiK 

C4H6,1butyne 

■Mi 

CaO+ 

1M\ 



Mn 


S03 

■p 

CH3 


C4H6,2butyne 


Ca(OH)2 


CuO 

Btir/il 





CH3C(CH3)2CH3 


C4H6, butadiene 


i«gt»i:iww=a-' 


CuO(cr) 






CH2CO,ketene 


C4H6,cyclo- 

168 

Ca(OH)2(L) 

224 

CuOH 

Bi?il 





Note: No subscripts or no phase for gaseous species; (cr): crystal; (a), (b), (c) (d), (1), (11): different crystal 
structures, (L): liquid, (gr): graphite. For detailed explanation of the species, see ref. 14. 


Free atom fraction is highly dependent upon experimental conditions and 
instrumentation. The common factors affecting the measured P are the flame chemistry, 
flame temperature and pressure, burner designs, nebulization efficiencies, observation 
height and entrainment of air in the flame which affects the fucl/oxidant ratio. Because of 
these reasons, investigators prefer to bypass these issues and prepare a calibration 
curve" , or report P values relative to a reference element (e.g. Cu, Mg, Ag) that may be 
assumed to be completely atomized (P = 1) at the flame condition. The review article 24 
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used in this report for comparison purposes, assumes that Cu is completely atomized in 
air-Cifh flame. 

bl. C 2 H 2 -0 2 -He system 

As mentioned before, the first case studied is the combustion of acetylene in oxygen 
diluted by helium. The calculated P values are reported in Table 4. 


Table-4 Reported vs. Calculated p Values, C 2 H 2 -Q 2 -He System 


O/F ratio = 1.368 



Calculated from 
CEA 2750 K 

Reported 24 
2750 K 

P Ratio 
calc ./reported 

Ca 

.573 

.22 

2.60 

Cr 

.763 

.41 

1.86 

Cu 

.999 

.84 

1.19 

Fe 

.985 

.82 

1.2 

O/F ratio = 1.758 



Calculated from 
CEA 2750 K 

Reported 24 
2750 K 

Ratio 

calc ./reported 

Ca 

.151 

.15 

E01 

Cr 

.245 

.25 

.980 

Cu 

.997 

1.04 

.959 

Fe 

.908 

.97 

.936 


Excellent agreement is seen between the calculated and predicted P values except for Ca 
and Cr at O/F ratio of 1.368. Even for this case, the ratio is well within a factor of 3.0. It 
is reasonable to say that the P values predicted from CEA would be safe to extend to the 
CH4 and RP-1 fuels for Ca, Cr, Cu and Fe. 

b2. H 2 -0 2 System 

Tejwani et. al. reported an occasion' when the material loss from bearing balls on the low 
pressure oxidizer turbo pump shaft in the space shuttle main engine could be measured 


Table-5 Reported vs. Calculated p Values H 2 -0 2 System 


O/F ratio = 6.0, pressure 2.455 atm. 


Calculated from 
CEA 3097 K 

Reported 2 
3097 K 

P Ratio 
calc./reported 

Cr 

.098 

.056 

1.75 

Cu 

.985 

.77 

1.28 

Fe 

.752 

.49 

1.53 

Mn 

.993 

.57 

1.74 

Ni 

.983 

.91 

1.08 
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after a flightworthiness test and P values could be estimated for some elements in the 
alloy. These values along with the values calculated from the CEA code are given in 
Table-5. 

In this case, Tejwani divided the plume cross section into three layers - one central layer 
and two outside layers. Only the central layer containing the Mach diamond and fairly 
constant temperature of 3097 K is considered here for comparison. 

Here the agreement is well within a factor of two. Considering the experimental error and 
many assumptions made in arriving at the final P values, this agreement is considered 
quite acceptable. 

c. Calculation of p values using CEA code 

Based on the results of Tables 4 and 5 above, it can be said that P values generated from 
the CEA codes for Ca, Cr, Cu, Fe, Mn, and Ni are fairly reliable. The value for Ag at the 
high temperature and fuel rich conditions should be fairly close to one - even though the 
CEA data base on Ag is not extensive. It is expected that P for A1 should also be reliable 
because the CEA data base on A1 species is fairly complete. The current data base for Co 
is not extensive so P for Co results should be treated with caution. 

Having thus generated some degree of confidence in the CEA data base and the free 
energy minimization method of determining free atom fractions, we now embark upon 
the prediction of P values for the nine elements under the chemical environment of direct 
interest to us. 

cl. CH4-O2 system 

The range of O/F ratio of 2.0 - 2.5 is used in this report because some of the initial 
experiments with the pilot thruster at present are being conducted at these ratios. It is 
likely that in future this range may be increased to cover O/F ratios from 1.75 to 4.0. New 
P values can be generated then. 

The data is also presented graphically in Figure- 1. The graphs are provided to show the 
extreme non-linear nature of P as a function of temperature for many elements. The 
temperature scale is the same in all graphs; however, the ordinate varies to accommodate 
the range of P values for a particular element. The ordinate scale for an element is kept 
the same for CH 4 -O 2 and RP-I-O 2 cases presented in the next section. 

c2. RP-I-O 2 with and without Sulfur 

RP-1 specifications allow for a maximum sulfur content of 30 ppm by mass. This case 
was run as a worse case scenario. It was observed that except for nickel, the calculated P 
values for all other metals did not change in the presence of sulfur. Even in the case of 
Ni, P values only at 1600 K were different - 1.000 without sulfur and 0.999 with sulfur. 
Co and Mn do not have compounds with S in the database; therefore, the effect of sulfur 
on their free atom fractions can not be stated conclusively. Apparently, the sulfur 
compounds of the remaining elements are not stable at the temperatures considered, thus 
they do not affect the free atom fractions. 
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Table-6 Calculated p Values for CH4-O2 System 


O/F ratio 2.0 

Temp K 

1600 

2000 

2400 

2600 

2800 

3000 

3100 

Ag 

1.000 

1.000 

1.000 

1.000 

0.999 

0.998 

0.995 

Al 

0.000 

0.000 

0.001 

0.003 

0.007 

0.016 

0.023 

Ca 

0.000 

0.002 

0.023 

0.053 

0.099 

0.149 

0.160 

Co 

0.559 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

Cr 

0.002 

0.458 

0.401 

0.380 

0.364 

0.356 

0.357 

Cu 

0.999 

0.998 

0.998 

0.997 

0.997 

0.996 

0.995 

Fe 

0.033 

0.490 

0.870 

0.911 

0.921 

0.918 

0.915 

Mn 

1.000 

1.000 

1.000 

1.000 

1.000 

0.999 

0.998 

Ni 

0.628 

0.999 

0.998 

0.998 

0.996 

0.995 

0.994 

Pd 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

O/F ratio 2.25 

Temp K 

1600 

2000 

2400 

2600 

2800 

3000 

3100 

Ag 

1.000 

1.000 

1.000 

1.000 

0.999 

0.998 

0.995 

Al 

0.000 

0.000 

0.001 

0.002 

0.005 

0.012 

0.017 

Ca 

0.000 

0.001 

0.013 

0.031 

0.062 

0.101 

0.118 

Co 

0.559 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

Cr 

0.001 

0.312 

0.269 

0.254 

0.244 

0.243 

0.249 

Cu 

0.998 

0.998 

0.997 

0.996 

0.995 

0.994 

0.993 

Fe 

0.019 



0.860 

0.880 

0.880 

0.879 

Mn 

1.000 



1.000 

1.000 

0.999 

0.998 

Ni 

0.628 



0.996 

0.995 

0.993 

0.992 

Pd 

1.000 



1.000 

1.000 

1.000 

1.000 

O/F ratio 2.50 

Temp K 

1600 

2000 

2400 

2600 

2800 

3000 

3100 

Ag 

1.000 

1.000 

1.000 

1.000 

0.999 

0.997 

0.995 

Al 

0.000 

0.000 

0.000 

0.001 

0.004 

0.009 

0.013 

Ca 

0.000 

0.001 

0.008 

0.019 

0.039 

0.069 

0.086 

Co 

0.559 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

Cr 

0.001 

0.195 

0.165 

0.156 

0.152 

0.158 

0.167 

Cu 

0.998 

0.997 

0.996 

0.995 

0.993 

0.992 

0.991 

Fe 

0.011 

0.237 

0.694 

0.791 

0.825 

0.833 

0.836 

Mn 

1.000 

1.000 

1.000 

1.000 

1.000 

0.999 

0.998 

Ni 

0.628 

0.999 

0.996 

0.994 

0.992 

0.990 

0.989 

Pd 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 
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Figure-1 P Values for CH4-02 System at Various Temperatures and O/F Ratios 
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Table- 7 Calculated p Values for RP-I-O 2 System 


O/F ratio 1.75 


1600 

2000 

2400 

2600 

2800 

3000 

3100 

Ag 

1.000 

1.000 

1.000 

1.000 

0.999 

0.993 

0.985 

Al 

0.000 

0.000 

0.001 

0.004 

0.009 

0.019 

0.028 

Ca 

0.000 

0.003 

0.033 

0.071 

0.114 

0.129 

0.117 

Co 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.999 

Cr 

0.017 

0.499 

0.429 

0.403 

0.383 

0.371 

0.368 

Cu 

0.999 

0.999 

0.998 

0.998 

0.997 

0.996 

0.995 

Fe 

0.056 

0.602 

0.903 

0.928 

0.930 

0.924 

0.920 

Mn 

1.000 

1.000 

1.000 

1.000 

0.999 

0.997 

0.994 

Ni 

0.999 

0.999 

0.998 

0.998 

0.997 

0.995 

0.994 

Pd 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

O/F ratio 2.0 


1600 

2000 

2400 

2600 

2800 

3000 

3100 

Ag 

1.000 

1.000 

1.000 

1.000 

0.998 

0.992 

0.984 

Al 

0.000 

0.000 

0.001 

0.002 

0.006 

0.012 

0.018 


0.000 

0.001 

0.014 



0.078 



1.000 

1.000 

1.000 

1.000 


1.000 



0.008 

0.290 

0.238 



0.210 


Cu 

0.999 

0.998 

0.997 

0.996 



0.992 


0.025 

0.390 

0.806 

0.861 

0.873 

0.870 


Mn 

1.000 

1.000 

1.000 

1.000 

0.999 

0.997 

0.993 

Ni 

0.999 

0.999 

0.997 

0.996 

0.994 

0.992 

0.990 

Pd 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

O/F ratio 2.25 


1600 

2000 

2400 

2600 

2800 

3000 

3100 

Ag 

1.000 

1.000 

1.000 

1.000 

0.998 

0.992 

0.983 

Al 

0.000 

0.000 

0.000 

0.001 

0.004 

0.009 

0.013 

Ca 

0.000 

0.001 

0.007 

0.017 

0.033 

0.049 

0.053 

Co 

1.000 

1.000 

1.000 

1.000 

1.000 

0.999 

0.999 

Cr 

0.004 

0.146 

0.115 

0.106 

0.103 

0.111 

0.121 

Cu 

0.998 

0.997 

0.995 

0.994 

0.992 

0.990 

0.988 

Fe 

0.013 

0.239 

0.675 

0.765 

0.792 

0.799 

0.803 

Mn 

1.000 

1.000 

1.000 

1.000 

0.999 

0.996 

0.993 

Ni 

0.998 

0.998 

0.995 

0.993 

0.990 

0.986 

0.985 

Pd 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 


The tabular data is presented in the graphical form in Figure 2. Once again, the non-linear 
nature of P is exemplified. 
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Figure-2 p Values for RP-1-02 System at Various Temperatures and O/F Ratios 
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d. P values for Pd 

As a trial case, P values were calculated with eight elements without Ag in the database. 
No difference in P values was seen for these elements compared to the case that included 
Ag. Apparently, C, H and O species are abundant enough in the mixture so that presence 
or absence of these elements in part per million quantities does not affect the result. 
Similarly, it is reasonable to assume that the absence of Pd from the CEA data base will 
not affect the P value for the nine elements documented in this report. 

Table-8 First Ionization Potential and Bond Energies (eV) of the Ten Elements 


At. No. 


I.P. 

M-O 

M-H 

M-S 

M-M 

13 

Al 

5.986 

5.31 

2.95 

3.88 

1.93 

20 

Ca 

6.113 

4.81 

1.74 

3.25 

0.16 

24 

Cr 

6.766 

4.43 

2.9 

4.14 

1.61 

25 

Mn 

7.435 

4.17 

2.4* 

3.12 

0.44 

47 

Ag 

7.576 

2.21 

2.34 

2.25* 

1.68 

28 

Ni 

7.635 

4.06 

3 

3.73 

2.71 

29 

Cu 

7.726 

3.55 

2.9 

2.95 

2.09 

27 

Co 

7.86 

3.81 

2.34* 

3.55 

1.73 

26 

Fe 

7.87 

4.24 

1.86* 

3.34* 

1.04 

46 

Pd 

8.34 

2.0 

2.43* 


0.00 


Values are taken from Lange’s Handbook 31 , except those marked with asterisk are from CRC Handbook 3 


A P = 1 is suggested for Pd at the temperature range under consideration. First ionization 
potential (I.P.) and bond energies for M-O, M-OH, and M-H bonds are given in Table 8. 
It can be seen that Pd possess the highest first ionization potential of all of the ten 
elements under consideration, while at the same time, Pd-O, Pd-H and Pd-OH bond 
energies are lowest of all M-O, M-H and M-OH bonds. Pd does not exist as a dimer 
(Pd 2 ). These facts indicate that the most likely stable specie at high temperatures would 
be Pd in its free atomic form (P = 1). 


Table-9 Estimates of Standard Heats of Formation of Pd Species 



MW 

H°f kcal/mole 

Pd (cr) 

106.42 

0 

Pd gas(30) 

106.42 

377 

Pd+ 

106.41945 

804.687 

Pd 2 H 

213.84794 

468.92 

PdH 2 

108.436 

468.92 

PdO 

122.4194 

192.971 


Therefore, it is recommended that a P value of one be used for Pd, until such time that 
full thermodynamic data become available on compounds of Pd. A few useful 
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references 27 ' 30 that contain infonnation on Pd species relevant to the hydrocarbon 
environment under consideration are included in the list of reference. 

If one were to neglect the entropy contribution, the heats of formation at 298.15 K for the 
following Pd species Pd(cr), Pd, Pd+, Pd 2 H, PdH 2 and PdO can be calculated from the 
bond energies reported in Table 8. These values are given in Table 9. These data may 
also be added to the CEA data base to provide a first order approximation of species 
distribution and an estimate of P, even though the temperature dependence data 
(coefficients ai - a 7 , bi, b 2 in equations 2 -4) on these species is not available. 

e. Addendum needed to CEA Data Base 

The following data on Co and Pd should be added to either the original CEA database or 
to the data base used by an individual: 

1 . Thermodynamic data on CoO (cr), C 03 O 4 (cr), and C 0 SO 4 (cr) as functions of 
temperature are available in the JANAF Tables. The nine thermodynamic 
constants used in equations 2, 3 and 4 will have to derived from this data before it 
can be inputted into the data base. Dr. Russell Claus at the Glenn Research 
Center has agreed to generate these coefficients. P values for Co may be updated 
once these species have been inserted in the CEA data base. Several other 
compounds of Co of interest to us (CoS 2 , C 03 S 4 , CogSg) are reported in the 
literature 30 . However, no useful thennodynamic data is available on these 
compounds. 

2. Add NiH and NEO3 to the data base. 

3. Add Pd species Pd (cr), Pd, Pd+, Pd 2 H, PdH 2 and PdO to the data base. 

4. Compounds of Ag of interest to us, such as AgOH, Ag 2 0, Ag 2 0 2 , Ag? 03 , Ag 3 04 , 
Ag 3 0, Ag 2 S reported in literature should also be added to the data base. It is 
suspected that many of these compounds may not be thermally stable at the high 
temperatures of interest to us, because the bond energies of Ag-0 and Ag-S are 
the lowest in the nine elements, and thus may not affect P for Ag. 

The CEA data base in its current form does not allow addition of a new element such as 
Pd and species based on it. 

Conclusion 

Reasonable conformity is seen between the P values reported in the literature and those 
calculated from the CEA program. P for six elements (Al, Ca, Cr, Cu, Fe and Ni) out of 
the ten elements may be used until actual experimental data become available. P = 1 for 
Pd is suggested based on the bond energy considerations. The calculated P values for the 
remaining three elements (Ag, Co and Mn) are suspected to be on the high side and may 
not be used with high degree of confidence. 
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Results from pilot tests with RP-1 on the CR-101 thruster can be used to establish the 
base line for the trace amounts of these elements in RP-1 fuel. The thruster does not have 
any moving parts, therefore, it can also be used for plume seeding experiments with these 
elements to generate P values and compare them with those predicted from the CEA 
code. 

The CEA data base should be extended with appropriate oxides, hydrides, hydroxide and 
sulfur containing species of Pd, Ag, Co and Mn as relevant thermodynamic data become 
available. 
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